Coupling-dependent rates of energy transfers from photoexcited Mott insulators to 

lattice vibrations 



Kenji Yonemitsu^^ and Nobuya Maeshima'^ 



0^ 
O 
O 
(N 

D 



u 

a 

o 
o 



00 

cn 
m 

(N 
O 
On 
O 



' Institute for Molecular Science, Okazaki 444-8585, Japan 
^ Department of Functional Molecular Science, Graduate University for Advanced Studies, Okazaki 444-S585, Japan 
Institute of Materials Science, University of Tsukuba, Tsukuba 305-8573, Japan 
(Dated: February 19, 2009) 

Photoexcited states are relaxed by transferring energy to the environments. In order to study 
which coupling allows fast energy transfer to lattice vibrations in correlated electron systems, we 
calculate the time evolutions of the kinetic energies of different types and frequencies of lattice 
vibrations. The one-dimensional half-filled Hubbard model is augmented with electron-lattice cou- 
plings that modulate transfer integrals, site energies, and Coulomb repulsion strengths. The time- 
dependent Schrodinger equation is solved for exact many-electron wave functions, and the classical 
equation of motion for the lattice displacements. In order to transfer energy to classical lattice 
vibrations that modulate transfer integrals or site energies, the translational invariance must be 
broken to give optical activity to an electronic excitation with wave number tt and to these lattice 
vibrations. On the other hand, a certain amount of energy is always transferred to lattice vibra- 
tions that modulate Coulomb repulsion strengths, irrespective of the symmetry of the ground state, 
as long as the corresponding electron-lattice couplings are present. In strongly correlated electron 
systems, these couplings can be strong, although they are usually insignificant because their effects 
on the equilibrium properties can be absorbed into redefinition of Coulomb repulsion strengths. We 
will discuss competition or collaboration between energy transfer pathways through different types 
of electron-lattice couplings. 



I. INTRODUCTION 

Photoinduced phase transitions have attracted much 
attention as intriguing cooperative phenomena in 
nonequilibrium conditions. ^'^ The electronic properties 
caused by many electrons are macroscopically changed. 
Especially in strongly correlated electron systems, they 
are changed often on an ultrafast time scale. The den- 
sity of photons needed for such a change is usually much 
lower than that of electrons involved because the energy 
supplied by photoirradiation is much higher than thermal 
energies. 

Quasi-one-dimcnsional Mott-Hubbard insulators 
such as halogen-bridged nickel-chain compounds'^ and 
Sr2Cu03^i^ have been expected as nonlinear optical 
materials with large nonlinear coefficients and fast 
response times. Ultrafast recovery times are desirable 
for optical high-speed switching devices. Then, fast 
energy dissipation from the photoexcited state into 
the environments is required. For Sr2Cu03, the ob- 
served photoinduced absorption is to an even-parity 
two-photon state that occurs immediately above the 
absorption edge, which is theoretically explained on the 
basis of the one-dimensional extended Hubbard model 
with alternating site energies.'* In such one-dimensional 
correlated electron systems, the spin-charge separation 
plays a crucial role in enhancing the nonlinear optical 
response 1^ 

Halogen-bridged nickel-chain compounds are known to 
exhibit a photoinduced insulator-to-metal transition and 
ultrafast decay of the photoinduced state. ^ Their pho- 
toinduced metallic properties are theoretically studied in 
the one-dimensional Hubbard modelj^ The difference be- 



tween the photoinduced state in the nickel-chain com- 
pounds and that in the copper-oxide chains is experi- 
mentally clarified by comparisons between their optical 
conductivity spectra and photoconductivity spectrai^ 

Quasi-one-dimensional Mott-Hubbard insulators 
are often accompanied by the dimerization that al- 
ternates transfer integrals to produce a spin-Peierls 
phase. Photoirradiation of such insulators as K- 
tetracyanoquinodimethane (K-TCNQ) are known to 
exhibit melting of the dimerization, which is called 
a photoinduced inverse spin-Peierls transition.^ Ul- 
trafast photoinduced melting of the spin-Peierls 
phase is observed with a very short decay time of 
photocarriersJ^ In bis(ethylenedithio)tetrathiafulvalene- 
difluorotetracyanoquinodimethane (ET-F2TCNQ) , 
which is regular without dimerization, photocarrier 
doping is concentrated on a Drude component and a 
midgap state is never formed owing to the spin-charge 
separationiii For a larger photoexcitation, the decay 
time of the metallic state is significantly shorter, suggest- 
ing that electron-electron scattering plays an important 
role. 

Recently, ultrafast charge dynamics in K-TCNQ, Rb- 
TCNQ, and ET-F2TCNQ with different magnitudes of 
electron-lattice interactions are compared In K-TCNQ 
and Rb-TCNQ, photocarriers are localized as polarons 
with about 70 fs and recombine with a few ps. In ET- 
F2TCNQ, photocarriers lead to a metallic state and de- 
cay with about 200 fs. These differences in the charge 
dynamics are reasonably interpreted as due to different 
magnitudes of electron- lattice interactions. The pho- 
toinduced inverse spin-Peierls transition is also theoreti- 
cally studied in the one-dimensional dimerized Hubbard 
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modeli^ and in its extension to include nearest-neighbor 
repulsionJ^ To understand photoinduced midgap states, 
however, electron-lattice interactions and the rapid for- 
mation of polaronic states are required for K-TCNQii^ 

In correlated electron systems on the one-dimensional 
regular lattice at half filling, the spin-charge separation 
holds in the low-energy limit and almost holds even in 
photocxcited states. Consequently, electron-lattice inter- 
actions play an important role in the decay of photocar- 
riers. In studying the relaxation of photoexcited states 
in correlated electron systems, it is numerically hard to 
observe the time evolution of the conductivity spectrum 
or the Drude weight, but it is easy to observe the time- 
dependent kinetic energies of phonons. In this paper, 
we then focus on the energy transfers from photoexcited 
electron systems to phonons. 

Quite recently, enhancement of phonon effects on an 
exciton has been reported in photoexcited states of one- 
dimensional Mott insulators on the basis of a density 
matrix renormalization group calculation for the one- 
dimensional extended Hubbard- Holstein model. Here, 
quantum phonons are incorporated and play an impor- 
tant role because the holon and the doublon are dressed 
with quantum phonons. For photoinduced dynamics near 
a quantum critical point, where quantum phonon fluctu- 
ations are essential, coherence is shown to be enhanced 
within the quantum Blume-Emery-GrifBths modelJ^ In 
this paper, however, we will limit ourselves to classi- 
cal phonons, which can suggest that electron-lattice cou- 
plings that are usually ignored can lead to ultrafast decay 
of the photoinduced state in correlated electron systems. 

We here briefly summarize the relation between quan- 
tum and classical phonons. Phonons originally behave 
quantum-mechanically. When the average number of 
phonons at each site is much larger than the fluctuating 
component, they behave classically. If phonon energies 
are not negligibly small compared with electronic kinetic 
energies, phonons follow the motion of electrons to pro- 
duce slowly moving polarons. The quantum nature of 
phonons must be retained to produce polarons in a trans- 
lation invariant manner. With decreasing phonon ener- 
gies, the motion of phonons is more retarded than that of 
electrons, and finally electrons feel as if phonons produce 
a quasi-static field. Phonons barely follow the motion of 
electrons to create hardly moving polarons. In the clas- 
sical picture, polarons break the translational symmetry, 
but quantum fluctuations restore the symmetry by the 
linear combination of differently placed polarons. In re- 
ality, thermal fluctuations would help the localization of 
polarons. 

In K-TCNQ, photocarriers are regarded as converted 
into polarons^"' presumably by relaxing the displace- 
ments that modulate transfer integrals fi-^ It is interpreted 
that photocarriers are localized and their recombina- 
tion (i.e., the decay of photocarriers) is slowed down by 
the corresponding electron-lattice interaction^^ The fact 
that the materials with stronger electron-lattice interac- 
tions have the longer decay times suggests either that 



electron-electron interactions are mainly responsible for 
the decay processi^ or that a new type of electron-lattice 
interactions compete with those responsible for the for- 
mation of polarons in determining the decay time. 

In this paper, we will suggest the latter possibility by 
introducing electron-lattice interactions that modulate 
Coulomb repulsion strengths. These interactions allow 
energy transfers to phonons, irrespective of the symme- 
try of the ground state. As electron-lattice interactions 
that modulate transfer integrals increase, the dimeriza- 
tion is enhanced, and the energy transfers to phonons 
that modulate Coulomb repulsion strengths are shown to 
decrease in this paper using classical phonons. This fact 
will survive quantum fluctuations of phonons. Numerical 
calculations for time evolution with quantum phonons on 
a reasonably large lattice at half filling are unfeasible and 
beyond the scope of this paper. 



II. MODEL AND METHOD 

In order to study the energy dissipation from a cor- 
related electron system into the environments, we take 
a one-dimensional half-filled Peierls-Holstein-Hubbard 
model with another type of electron-lattice couplings that 
modulate on-site repulsion strengths, 
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where c j (cj,^) creates (annihilates) an electron with 



spin a at site j, Uj^a = Cj^^cj^a, nj = J2cr^j,<r> *o de- 
notes the bare transfer integral, and U the bare on-site 
repulsion strength. As for electron-lattice couplings, we 
consider different types and frequencies of lattice vibra- 
tions. The lattice displacements "j"'^ ^^'^ ^j™^ 
modulate the transfer integral, the site energy, and the 
on-site repulsion strength, with coupling strengths, sL'"'', 
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and s-y , respectively, at bond or site j. Here, s 



(m) 



(i^=a, (3, 7) are not smaller than zero, and the suffix m 
denotes modes having the bare phonon frequencies, to, 
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The quantities ii 
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, and Wj are 

the time derivatives of u""^^ , ^'j™'* i ^j™"* ; respectively. For 
simplicity, these displacements are regarded as indepen- 
dent. 
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We define each type of electron-phonon coupling 



strength by Aq 



and 



As long as the lattice vibrations 



are treated classically, the ground state is given by 
=0. In this case, the ground state is 



u 



(rn) ■ (m) 



determined not by the distribution of {sq™''}, {sJj'"''}, or 



{s^^™''}, but by their sums, Aq, A^, and A^, as derived in 
a straightforward manner from the Hellmann-Feynman 
theorem. Therefore, we will specify the electron-phonon 
coupling strengths by Aq,, A^j, and A-y in the following, al- 
though the dynamics induced by any excitation depends 
on the distribution of couplings. 

In this paper, we consider many phonon modes of dif- 
ferent frequencies for each type of electron-phonon cou- 
plings to reduce coherent oscillations of energy exchange 
between the electronic and lattice systems. Therefore, we 
take a wide distribution of phonon frequencies: tJ\^^ — 
iOi, X m/100 for v=a^ /?, 7 and m—1, ■ ■ •, 100 with uj^ be- 
ing the maximum phonon energy uj^, = max^ w^.™"* . Fur- 
thermore, we take a uniform distribution of couplings: 
s^J"'* — The number 100 is more than realistic val- 
ues even for complex molecular materials, but it is so 
taken as to be large enough in that the numerical results 
are almost unchanged if we take 10* instead. Although 
phonons are treated classically in this paper, we employ 
exact many-electron wave functions for the ground state 
and the photoinduccd dynamics. 

Photoexcitations arc introduced by adding 
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X < exp 



lea 



dtE{t) 



- 1 ct c,+i,, + H.c. (2) 



to Eq. ( [T]) , where e is the absolute value of the electronic 
charge, a is the lattice spacing, and e, a, and Ti are set 
to be unity. The time-dependent electric field E{t) is 
given by E(t) — E'oxt sinwoxt^, with amplitude i?oxt and 
frequency Woxt for Q < t < T^^ [E{t) is zero otherwise.] 
with Tirr = 27rA'^oxt/woxt and integer iVoxt- 

The time-dependent Schrodinger equation for the ex- 
act many-electron wave function | V-'(^)) is numerically 
solved by 



^]{t + dt)) ~ exp 
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where the time evolution operator is expanded as 
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with time slice dt— 10 ^ to the (n=)15th order and by 
checking the conservation of the norm and of the total 
energ y f or t > Ti„. This method is the same as in 
Refs. Hllil but we here use a much smaller value for 



dt to raise the precision. The classical equation of mo- 
tion for the lattice displacements is solved by the leapfrog 
method, where the force is derived from the Hellmann- 
Feynman theorem. 

To observe the energy dissipation into the environ- 
ments, we calculate the time evolution of the kinetic en- 
ergy of each type of phonons. 
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and compare the time evolution of the total energy, 
Etot = (H). In the following, the distortion amplitudes 
in the initial ground state are defined by the dimension- 
less quantities. 
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III. RESULTS 

In this paper, we use to=l and U=& for Mott insula- 
tors and cJcxt^S, which is slightly above the lowest optical 
excitation energy. We take A^ext=20, so that riir=41.9, 
which corresponds to 27.5 fs (275 fs) if to is 1 eV (0.1 eV). 
Furthermore, we use £'ext=l, which is not so large in that 
the increment in the total energy ranges about 1.5 to 3 
(i.e., 0.5 to 1 photon) in the present 12-site system. Al- 
though it is not easy to obtain general results on the dif- 
ferences between Mott and band insulators, we first com- 
pared the energy transfer to phonons in a Mott insulator 
and that in a band insulator, where the bare transfer in- 
tegral tQ in Eq. ([T]) is replaced by t^ — {—lySt, 5t is set at 
0.75, and [/ at to have a charge gap of 3. For A^=l and 



= 1, the ratio of the increment in £'kin({w}) 



to that in E'tot, Aii^kin({w})/A£'tot, increases to about 
8x10-2 and A£'ki„({u})/A£;tot to about SxlO'^ for the 
distortion amplitude uqs to 0.5. They are one order of 
magnitude larger than AEkin{{u}) / AEtot < 3x10"'^ and 
A£'kin({w})/A£'tot < 8 X lO"'' in a Mott insulator with 
dimerization (i.e., a spin-Peierls insulator) even with a 
larger A^ and the same phonon frequencies. 
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We have calculated same quantities for another type 
of band insulators where the site energies are alternated 
instead of the transfer integrals. In such a case also, the 
energy transfer to phonons is generally larger than that 
in a Mott insulator. This can be understood as follows. 
In principle, non-interacting electrons in the translation 
invariant system cannot have a charge gap at half filling. 
In order for them to have a gap, they must be accompa- 
nied by alternation of transfer integrals, site energies, or 
other potentials, either intrinsically by electron-phonon 
couplings or extrinsically by St, for instance. Then, the 
originally inactive (i.e., in the translation invariant sys- 
tem) electronic excitation with wave number tt gets opti- 
cal activity. It means that a certain amount of energy is 
always transferred to phonons that are coupled with this 
excitation. Meanwhile, a Mott insulator has a gap even 
in the translation invariant system, where this excitation 
remains optically inactive. This suggests that another 
type of clcctron-phonon coupling is effective for a Mott 
insulator to allow large energy transfer to phonons, com- 
pared to a band insulator with a similar magnitude of 
gap. 

From now on, we focus on Mott insulators. Although 
we use a particular set of parameter values to show nu- 
merical results, qualitative aspects are general and ob- 
tained for other parameter values as well. As an example, 
Figs.[Ha) and[TJb) show the time evolution of -Ekin({u}) 
and that of E'kin ({w}) , respectively. The bare phonon fre- 
quencies adopted here range from 0.01 to 1, so that their 
periods range from 6.3 to 628. Thus, the rapid oscilla- 
tions (of period about 2) are due to either forced or al- 
most resonantly-induced electronic excitations, which are 
coupled with these phonons. Some structures on larger 
time scales [of periods about 6 and 12 in Figs. [T](a) and 
[Ijb)] are also visible. Their appearance is allowed by 
the couplings with phonons. In the figures, the arrows 
indicate the pulse duration, Tjir. After the oscillating 
electric field is switched off, t > Tirr, the forced oscilla- 
tion disappears, but rapid oscillations due to electronic 
excitations survive. The energy transfer to phonons con- 
tinues on a large time scale [of the largest period 628 
(not shown)]. It becomes conspicuous especially for low 
(maximum) phonon frequencies. The dependence on the 
phonon frequencies is discussed later. 

For Aq values used in these figures, the ground state 
is in the spin-Peierls phase that has alternating u^™^ oc 
(—1)-'. Without this dimerization, the ground state is 
translation invariant and E]^[n{{u})—Ekin{{v})—0, which 
is numerically confirmed in the present finite-size sys- 
tem with small Aq.. It should be noted here that this 
is not really the case if phonons are treated quantum- 
mechanically as in Ref. \W, where the holon and the 
doublon are dressed with quantum phonons. As long 
as phonons are treated classically, the wave function is 
a direct product composed of the electronic and lattice 
parts where electrons feel the averaged phonon density 
and phonons feel the averaged electron density so that 
electron-phonon correlation is missing. Then, the photo- 
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FIG. 1: (Color online) Time evolutions of kinetic energies 
of (a) phonons that modulate transfer integrals, -E'kin({ii}), 
and (b) phonons that modulate on-site repulsion strengths, 
Eii_in{{w}), with different coupling strengths Xa- The other 
parameters are to=l, U=6, A/3=3, X'y=2, u>ct=i^i3~u!.y — l, 
idc^t~S, Ncxt~20, and i5ext=l. The arrows indicate the pulse 
duration, Tin- 



carriers cannot be dressed. Nevertheless, we expect that 
the different characters of the electron-phonon couplings 
considered here basically survive quantum fluctuations. 
Unless X.y=0, £'kin({w}) always becomes finite because 
the optical excitation involves the modulation of double 
occupancy {n-j.-^rij^i), irrespective of whether the ground 
state is dimerized or regular. Figure [TJa) shows that 
E]iini{u}) increases with Aq. It is because Aq increases 
the distortion amplitude uqs , which enhances the optical 
activity of the electronic excitation with wave number tt 
by deviating further from the translation invariance. Fig- 
urc[T]Jb) shows that E]^in{{w}) decreases as Aq increases. 
It is because Aq and thus uqs increase the antiferromag- 
netic spin correlation on the strong bonds (with the larger 
transfer integrals), which makes the modulation of dou- 
ble occupancy harder. These facts will be made clearer 
when they are plotted as a function of the distortion am- 
plitude later. 

Not all the energy transfer to phonons is regarded as 
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FIG. 2; (Color online) (a) Ratio of increment in -Bkin ({""}) to 
that in iJtot, as a function of the distortion amplitude uqs, 
with different coupling strengths A-y . (b) Ratio of increment in 
i'kin({w}) to that in -Btot, as a function of the distortion am- 
plitude wgs, with different coupling strengths Aa. Both are 
measured at t=100. The parameters are to=l, U=6, A^=3, 
LOa—LLii3=iO'j=l, ajcxt=3, A''cxt=20, and iJcxt = l- 



FIG. 3: (Color online) Ratios of (a) increment in i5kin({^'}) to 
that in i?tot, and (b) increment in i?kin({v}) to that in -Etot, 
as a function of the distortion amplitude uqs , with different 
maximum phonon frequencies ijj=ujct—uJis=ujry. Both are mea- 
sured at t=100. The parameters are to — l, U—6, A0=3, A^=2, 
Wext=3, Afext=20, and £ext=l. 



contributing to dissipation. Some of the transferred en- 
ergy will flow back into the electronic system and forth 
into the lattice system, causing oscillations of kinetic 
energies of phonons. To estimate the component that 
contributes to dissipation, in other words, to remove 
the effect of such energy oscillations, we draw a curve, 
for instance for _Ekin({'i*}), by connecting local minima 
in time (when it oscillates), measure its value [which 
is not greater than the bare £'kin ({«})] at the observa- 
tion time, and regard it as the increment in E\.i^{{u}), 
A£'kin({u}). Figure mja) shows the ratio of A£'kin({u}) 
to the increment in the total energy AE'tot , as a function 
of the distortion amplitude ugs- In the classical pic- 
ture of phonons, unless the ground state has a finite lat- 
tice distortion, the lattice remains undistorted even after 
the photoexcitation, which does not break the transla- 
tional symmetry. That is, for itGS=0, no energy is trans- 
ferred to w-phonons. As a consequence, uqs is a control- 
ling parameter for the energy transfer to u-phonons. As 



noted previously, it is a monotonically increasing func- 
tion. This quantity increases with the coupling A-y also. 
The coupling X-y reduces the on-site repulsion strengths 
and enhances the superexchange interactions between the 
spins on the strong bonds, which strengthens the effective 
coupling between electrons and w-phonons. Figure (Hb) 
shows the ratio of Ai^kin({w}) to A£^tot, as a function of 
the distortion amplitude WQg. It is a steeply increasing 
function above a threshold. The threshold is present be- 
cause the distortion amplitude wgs defined in Eq. (fTU]) 
is finite in the A-y — > limit. This quantity is suppressed 
by the coupling Aq , as explained previously. Namely, the 
energy flow into w-phonons is suppressed if the dimeriza- 
tion of transfer integrals is large. 

We vary the distributions of phonon frequencies to 
see their effect on the energy transfer to phonons. 
Figure [3fa) shows that the ratio of AE\iin{{u}) to Ai^tot 
increases with shifting the distributions of phonon fre- 
quencies to the low side. Here, the energy ratios are ob- 
served at i=100 shortly after the oscillating electric field 
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is switched off. If we decrease the phonon frequencies fur- 
ther, even the lattice vibration of the highest frequency 
LOa has the period longer than the observation time, lead- 
ing to decreasing of the energy ratio. However, if we 
observe it much later so as to allow a few cycles of the 
fastest lattice vibration before the observation time, the 
energy ratio tends to increase further. Although we can- 
not reach the adiabatic limit of w ^ due to numerical 
limitations, the energy flow into u-phonons seems to be- 
come the largest in this limit. This behavior is in contrast 
to that of w-phonons. Figure El^b) shows that the ratio 
of Ai?kin({w}) to AiJtot is almost completely suppressed 
if phonon frequencies are low. Note the difference in the 
scales of the ordinate axes. This energy ratio increases 
with shifting the distributions of phonon frequencies to 
the high side as long as phonon energies are lower than 
electronic excitation energies. If some phonon energies 
become comparable with electronic excitation energies, 
this quantity would show a complex behavior as a func- 
tion of phonon frequencies. 

These contrastive behaviors of u- and u-phonons are 
understood in the following way. The dimerization 
uj™' cx (—1)-' breaks the translation invariance, giv- 
ing the optical activity to an electronic excitation with 
wave number tt. This electronic mode is photoexcited 
almost resonantly here and transfers energy directly to 
M-phonons. However, the energy transfer to w-phonons 
is indirect. This electronic excitation makes the charge 
densities at even and odd sites different, which applies a 
force to u-phonons. It then induces oscillations of charge 
densities and those of w-phonons. The energy flow into 
thus indirectly coupled u-phonons is facilitated by mak- 
ing different excitations closer in energy. Consequently, 
the energy flow into w-phonons seems to become largest 
when their energies are comparable with those of elec- 
tronic excitations. 



IV. SUMMARY 

The energy flow from a correlated electron system 
into the environments after photoexcitation is theoret- 
ically studied in a one-dimensional half-filled Peierls- 
Holstein-Hubbard model augmented with another type of 
electron-lattice couplings that modulate Coulomb repul- 
sion strengths. We consider different types and frequen- 
cies of lattice vibrations, which are treated classically. 
The time-dependent Schrodinger equation is numerically 
solved for the exact many-electron wave functions. The 
time evolutions of the kinetic energies of different types of 
phonons are observed and compared with the increment 
in the total energy. When Mott and band insulators with 
a similar magnitude of gap are compared in the present 
half-filled model, the energy transfer to phonons is gen- 
erally larger in a band insulator than in a Mott insulator. 
Because the translational symmetry is broken to have a 
finite gap in the band insulator at half filling, the opti- 
cal activity is given to an electronic excitation with wave 



number tt and phonons coupled with it. 

In the spin-Peierls phase, the dimerization that alter- 
nates transfer integrals gives the optical activity to the 
electronic excitation with wave number tt, which allows 
the energy transfers to phonons that modulate trans- 
fer integrals and to phonons that modulate site energies. 
These energy transfers increase with dimerization ampli- 
tude, but their dependences on phonon frequencies are 
quite different. The transfer to phonons that modulate 
transfer integrals increases with decreasing phonon fre- 
quencies, while that to phonons that modulate site ener- 
gies increases with increasing phonon frequencies as long 
as these phonon energies are lower than the electronic ex- 
citation energies. This difference originates from the fact 
that the coupling with the optically active electronic ex- 
citation is direct for the former phonons and indirect for 
the latter phonons. Nevertheless, these energy transfers 
are generally smaller than those in band insulators. 

The above facts suggest that another type of electron- 
phonon couplings that modulate Coulomb repulsion 
strengths are effective in transferring energy to phonons. 
It should be noted that Coulomb repulsion strengths 
are generally much larger than the charge gap. There- 
fore, this type of electron-phonon couplings can be quite 
strong in principle. They are not usually regarded as 
substantial because they do not change the symmetry of 
the electronic state. Their effects can be absorbed into 
redefinition of Coulomb repulsion strengths as long as the 
equilibrium properties are concerned. Their effects can 
be significant in such nonequilibrium conditions as in the 
photoinduced dynamics. The energy transfer to these 
phonons always takes place irrespective of whether the 
ground state is dimerized or regular. They are relatively 
important when the dimerization is weak. 

In view of the energy transfer to phonons, those which 
modulate transfer integrals compete with those which 
modulate on-site repulsion strengths. With increasing 
dimerization amplitude, the antiferromagnetic spin cor- 
relation is strengthened on the strong bonds, making the 
modulation of double occupancy harder. It is interest- 
ing that the electron-phonon couplings that modulate 
on-site repulsion strengths enhance the energy transfer 
to phonons that modulate transfer integrals by reducing 
the on-site repulsion and thus by increasing the superex- 
change interaction. Above all, electron-phonon couplings 
that modulate Coulomb repulsion strengths can gener- 
ally be important when the energy dissipation from the 
photoexcited state is considered. Although we do not 
study here, the nearest-neighbor repulsion would also be 
modulated by relevant phonons and can be quantitatively 
important. For instance, in transition metal oxides, the 
repulsion between an electron in a d-orbital and another 
in an O p-orbital would be strongly modulated by the 
displacement of the oxygen ion. Because the oxygen ion 
is light and it can oscillate with high frequency, it would 
allow quite a rapid decay of the photoexcited state. 
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